Accelerated expansion and the virial theorem 
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When dark matter structures form and equilibrate they have to release a significant amount of 
energy in order to obey the virial theorem. Since dark matter is believed to be unable to radiate, 
this implies that some of the accreted dark matter particles must be ejected with high velocities. 
These ejected particles may then later hit other cosmological structures and deposit their momentum 
within these structures. This induces a pressure between the cosmological structures which opposes 
the effect of gravity and may therefore mimic a cosmological constant. We estimate the magnitude 
of this effect and find that it may be as large as the observed accelerated expansion. Our estimate 
is accurate only within a few orders of magnitude. It is therefore important to make a much 
more careful calculation of this redshift dependent effect, before beginning to interpret the observed 
accelerated expansion as a time dependent generalization of a cosmological constant. 

PACS numbers: 98.62.Mw, 98.62.Py, 98.80.Es 



I. INTRODUCTION 

The expansion of the universe has been observed to 
accelerate. This was first noted through the analysis of 
supernovae P, 0] and more recently using cosmic mi- 
crowave background and baryonic acoustic oscillations 
Most analyses have demonstrated that a cosmo- 
logical constant is in good agreement with data [1, , 
however, one may easily imagine models which could al- 
low for generalizations beyond a simple cosmological con- 
stant. Before we start parametrizing generalized models, 
or fitting new free parameters, it is important to carefully 
consider and remove known effects. 

One effect, which has not received much attention in 
this context, arises from the merging between dark mat- 
ter structures and is always present in a bottom-up for- 
mation scenario of cosmological structures. This effect 
induces an effective pressure between cosmological struc- 
tures, which in principle may lead to systematic errors in 
the interpretation of the cause of the acceleration. We 
will now discuss the origin of this effect. 

When a dark matter structure is formed and equili- 
brated, either through merging or accretion, it has to 
obey the virial theorem. The virial theorem quantifies 
the connection between the potential energy and the ki- 
netic energy of the entire structure, 2K + W = Q, where 
K is the total kinetic energy, and W is the total po- 
tential energy For instance, if a structure is formed 
by bringing cold particles in from infinity (i.e. initially 
J^init = Winit = 0), then the virial theorem implies that 
an energy corresponding to the resulting total kinetic en- 
ergy, K, must be disposed of For non-radiative dark 
matter particles this means that a significant number of 
particles typically are ejected with high momentum. 

In this Letter we explain and estimate this effect. We 
find that it indeed may have a significant contribution 
to the redshift dependent acceleration of the expansion. 
Since this is an unavoidable effect, we emphasize the need 
for a careful calculation of its affect on the cosmological 
expansion. 



II. EFFECTIVE REPULSION 

Let us consider a system of 3 equilibrated dark matter 
structures, 2 large and 1 small. The 2 large structures are 
placed at a large distance from each other, for instance 
at 100 times their virial radius. The small structure, 
called A, is placed close to one of the large structures, 
called B, for instance at 5 times the virial radius. A is 
given the circular velocity such that it is orbiting B in a 
stable circular orbit. The total force on the other large 
structure, called C, is given by the total gravity from A 
and B. 

Now, let us consider a slightly different configuration, 
namely one where A is given the same speed as before, 
but directed towards B. This implies that after a few dy- 
namical times A is essentially engulfed by B, and the 2 
structures A and B will reach a new equilibrium config- 
uration. The particles originating from A now sit at a 
deeper potential, and they therefore have higher kinetic 
energy on average. 

At the end of this equilibration the virial theorem must 
hold, and therefore approximately half the total change 
in energy must have been disposed of. If the structures 
are composed of coUisionless and non-radiative particles, 
then that energy can only be radiated away by "sacrific- 
ing" some of the incoming particles, or by increasing its 
size slightly. The ejected particles must leave the system 
with velocities above the escape velocity. These parti- 
cles will generally leave the system in a wide cone along 
the axis of collision, however, for this discussion we just 
assume that they leave the system spherically. 

What happens with these sacrificed particles with pos- 
itive total energy? Some will just leave the system ra- 
dially, and be of no further concern for now, however, 
some of them will happen to be directed towards the 
other large structure, C, and when they hit C they will 
transfer their momentum to it through standard gravita- 
tional effects like dynamical friction. The effect on C is 
therefore an effective pressure. 

When we compare the 2 configurations above, then 
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we see that in both cases there will be a gravitational 
acceleration of C towards the combined system of A and 
B. However, in the second configuration there will be an 
additional pressure on C. If we knew nothing about the 
merger history of A and B, then we might interpret this 
extra acceleration of C as a negative gravitational effect. 
We will refer to this as the "rejected acceleration". 

If we were concerned with measuring the actual accel- 
eration of the expansion of the universe, then it would 
be important to consider the magnitude of the effect de- 
scribed above, and if the effect would be non-negligible, 
then a careful subtraction must be done before we can as- 
cribe the acceleration to e.g. a time varying cosmological 
constant. 



III. HOW LARGE IS THE EFFECT 

We now wish to make a rough estimate of the effect. 
Let us consider the universe today, and let us assume that 
half of the total mass has been assembled into structures, 
all with the mass of a galaxy, M = IO^^Mq. Naturally 
one should consider the full mass distribution, N{M,z), 
however the argument below gives the same result for 
structures of lO^M© and within a factor of a 

few (largely from the difference in mass-concentration), 
so a proper integral over N{M, z) is expected to give the 
same result within an order of magnitude. 

Using a critical density of pc = 1.4 x lO^^Mo/Mpc'^, 
this means that the average distance between 2 structures 
is about 2.2 Mpc. 

We wish to compare the rejected acceleration (that is, 
the effective acceleration induced by the sacrificed par- 
ticles) to the acceleration which is generally believed 
to originate from a cosmological constant. The cos- 
mological constant induces an acceleration of the order 
GM{^{r)/r'^ ss 9 x 10"^^ m/s^ on a distance of r = 2.2 
Mpc, where G is the gravitational constant, and Mf^{r) 
is 0.7 pc integrated over a sphere of radius 2.2 Mpc. This 
simple estimate is, within a factor of 2, the same as one 
obtains when numerically solving the Friedmann equa- 
tion. We therefore have that the acceleration from a cos- 
mological constant is given by 

aA(2.2Mpc) « 9 X 10~" m/s^ . (1) 

Next we consider the acceleration from the sacrificed 
particles. First it is important to remember that there 
is a time delay from the time when the particles were 
merged onto (and ejected from) one galaxy, and till they 
were absorbed by another galaxy. This retarded time de- 
pends on the typical distance between structures. If the 
ejected particles from a galaxy has of the order 200 km/s 
(similar to the peak dispersion in the galaxy) then a dis- 
tance of 2.2 Mpc imphes that the particles were ejected 
roughly 10^° yrs ago. We therefore see that only the 
particles ejected with at least 200 km/s (in addition to 
the kinetic energy to leave the potential) will reach the 
other galaxy today. The ejected energy may easily be 



even higher or lower for some of the ejected particles. 
Numerical tests of the ejection mechanism seems to indi- 
cate that of the order 20 — 30% of the particles are ejected 
in a spherical cold collapse simulation [ll|, and between 
10% and 40% in merger simulations (l3 |. 

Let us now discuss the rate at which the particles 
are ejected from the first galaxy. The merging rate has 
clearly not been a constant throughout the history of 
the universe, however, to get the order of magnitude 
we just simplify by a linear merging rate in time, such 
that 5M/5t = M/tn = IO^'^Mq/IS.I Gyr. We there- 
fore use SM/6t = lO^^Afo/Myr. To fulfiU the virial the- 
orem, half the changed energy will be ejected or used 
to increase the size of the system. The distribution is 
merger dependent, however, for this order of magnitude 
estimate, we simply consider a quarter of the energy to 
be ejected which implies that roughly 0.25 SM /St will be 
ejected with v = 200 km/s. The acceleration is therefore 
dv/dt = v/M X 0.25 5M/St w 1.8 x lO^^^ m/s^. 

We finally have to consider how many of the ejected 
particles will deposit their momentum in another galaxy. 
If we consider a very long timescale, then all particles 
will eventually hit some galaxy. Equivalently, the one 
receiving galaxy may absorb particles from a range of 
previously merging structures. As a theoretical upper 
limit we therefore have, that all the ejected particles may 
be absorbed. However, many of these galaxies will lie 
at even larger distances, and the particles will therefore 
deposit their energy far into the future. In that case 
the velocities will be redshifted to lower values, and less 
energy will be deposited. Typically, the velocity must 
be redshifted below the escape velocity of the receiving 
structure p^ . 

We therefore conclude that the rejected acceleration 
must lie below the value found above 

arcjcct(2.2Mpc) < 1.8 x 10"" m/s^ . (2) 



IV. DISCUSSIONS 

We have seen above that the rejected acceleration is 
smaller than 1.8 x 10~^'^m/s^, which happens to be of 
the same order of magnitude as the observed acceleration, 
flA ~ 9 X 10"^^ m/s^ within a factor of a few. 

We are here entertaining the view of the redshift being 
a true Doppler shift, since we are ascribing the accelera- 
tion to a changed velocity of the galaxies and other cos- 
mological structure. This is possibly completely equiv- 
alent to the view of the expansion being a stretching of 
space |l3|. This question is still actively debated (see a 
list of references in and we will not enter that discus- 
sion here. Similarly, it is also actively debated how the 
differences in an inhomogenneous space between a local 
accelerating space and a globally accelerating space will 
manifest themselves observationally (see e.g. (Tsl. Kd^). 

We estimated above that the simplification from using 
one mass only gives a order of magnitude uncertainty. In 
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principle this can be calculated more accurately, by using 
the correct mass distribution, N{M,z), e.g. from Press- 
Schechter theory. One complication that will arise is that 
particles ejected from high mass structures are difficult 
to absorb by smaller structures due to the higher particle 
velocities. In that sense there may be a larger effect of 
the rejected acceleration on larger structures. 

The merging rate is both mass and redshift dependent, 
and therefore the merging rate (and hence the rejected 
acceleration today) may be somewhat different than the 
estimate used above. 

Another point which was ignored above is that the 
mass of the structures were smaller at the time when the 
particles were ejected. For instance the lO^^M© structure 
considered above would only have about have 1/4 of that 
mass 10^° yrs ago. We have also been using that a quarter 
of the energy must be ejected, however, if the first galaxy 
changes its size more than assumed during the merging, 
then it may hold more energy, and fewer particles need 
to be ejected, and vice versa. 

A rather non-trivial point is the retarded time: the sac- 
rificed particles will only be absorbed in other structures 
at a later time. This retardation will be different for dif- 
ferent particles, depending on the mass of the structure 
they are first merged onto. The reason is that whereas 
small mass structures are formed earlier (hence giving 
a larger retarded time) then the sacrificed particles are 
ejected at lower velocities. This all implies that there 
should be a strong redshift dependence on the rejected 
acceleration. It is also interesting to note that the effect 
of rejected acceleration will diminish in the far future, 
when structure formation has been reduced for a long 
time, and all ejected particles have either been absorbed 
or their energies redshifted away. 

To a first approximation this effect should be isotropic, 
as the observed acceleration is. However, there may be 
a slight anisotropic pressure near a large over-density 
where most of the merging happens along the filament. 
This possibility can be quantified by numerical simula- 
tions, by observing to which degree the ejected particles 
are emitted spherically or along the axes of merging. 

One might ask if this effect should already have been 
observed by cosmological N-body simulations, and the 
answer is no. Most simulations are executed by separat- 
ing the expansion of the universe from the local effects of 
structure formation. This means that the rejected accel- 
eration cannot induce an extra acceleration of the expan- 
sion, but at best impose an extra radial pressure on all 
structures in present day simulations. This radial pres- 
sure is low, and will be very difficult to disentangle from 
the dynamics of the normally infalling particles. It is not 
known to me if a method to perform cosmological sim- 
ulations exists, where one can dynamically include the 
effect of the cosmological constant: if it will be possible 
to include for instance a local density dependent acceler- 
ation due to a generalized cosmological constant, then it 



will certainly also be possible to include the effect of the 
rejected acceleration, at least in principle. 

Finally, a different point related to the retarded time is 
the velocity with which the particles are ejected from the 
original structures. There will naturally be a broad dis- 
tribution of velocities, since not all particles are ejected 
at the same velocity. One therefore has to fold this ve- 
locity distribution of the ejected particles with the mass 
distribution, to find the full redshift dependence of this 
acceleration. One should also include the effect that dif- 
ferent structures emit particles at different time, which 
may lead to an accumulation of the effect. For instance a 
IO^^Mq structure will reject particles earlier and slower, 
than a 2 X lO^^M© structure. Therefore, there may be an 
overlap of the arriving particles today, from the particles 
ejected at different times. Furthermore, above we merely 
used the proper distance today, however, in reality the 
relevant distance will be slightly shorter because of the 
expansion of the universe. Considering all these effects 
is a somewhat non-trivial calculation which we intend to 
address in the near future. 



V. CONCLUSION 

The observed accelerated expansion of the universe is 
generally believed to be caused by a cosmological con- 
stant, or some time dependent generalization hereof. In 
this brief note we point out a simply dynamical mech- 
anism which may have a significant contribution to the 
acceleration. 

The dynamical effect arises from merging of dark mat- 
ter structures, where the virial theorem implies that a 
significant amount of kinetic energy must be carried away 
by ejected particles. These ejected particles will later de- 
posit their momentum on other cosmological structures, 
inducing an effective accelerated expansion. 

We estimate the effective accelerated expansion from 
the ejected particles to be of the same order of magni- 
tude as the observed accelerated expansion of the uni- 
verse. The simple estimate presented here is accurate 
only within a few orders of magnitude. 

We emphasize that this dynamical effect is a well 
known physical effect which cannot be avoided. It is 
therefore something which needs to be calculated accu- 
rately and accounted for before measurements of a time 
varying cosmological constant will be trustworthy. 
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